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Intercombination Effects in Resonant Energy Transfer 
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We investigate the effect of intercombination transitions in excitation hopping processes such as 
those found in Forster resonance energy transfer. Taking strontium Rydberg states as our model 
system, the breakdown of LS'-coupling leads to weakly allowed transitions between Rydberg states 
of different spin quantum number. We show that the long-range interactions between two Rydberg 
atoms can be affected by these weakly allowed spin transitions, and the effect is greatest when there 
is a near-degeneracy between the initial state and a state with a different spin quantum number. 
We also consider a case of four atoms in a spin chain, and show that a spin impurity can resonantly 
hop along the chain. By engineering the many-body energy levels of the spin-chain, the breakdown 
of LS coupling due to inter-electronic effects in individual atoms can be mapped onto a spatial 
separation of the total spin and the total orbital angular momentum along the spin chain. 

PACS numbers: 32.80.Ee, 34.20.Cf, 37.10.Jk, 87.15.Hj 


I. INTRODUCTION 

Non-radiative exchanges of energy mediated by dipole- 
dipole interactions play a crucial role in a variety of pro¬ 
cesses, ranging from photosynthesis in natural biological 
systems [X; to highly efficient light emission in organic de¬ 
vices [2] . The usual requirements for a significant energy 
transfer are an electric dipole-dipole interaction between 
donor and acceptor molecules, and a near degeneracy be¬ 
tween the initial and final states that ensures the process 
is always resonant [3] . Resonant energy transfer has also 
been extensively studied in atomic physics, where these 
conditions are easily met. For example, dipole-dipole in¬ 
teractions in dense, optically excited samples can lead to 
cooperative Lamb shifts By using Rydberg states, 

rather than low-lying electronic states, the strength of 
the interaction can be increased by many orders of mag¬ 
nitude, leading to energy exchange over macroscopic dis¬ 
tances 

It is often assumed, in agreement with the electric 
dipole selection rules, that only states with the same 
value of the total electron spin quantum number S are 
coupled through dipole-dipole interactions. This assump¬ 
tion is justified for interactions between alkali atoms, 
since S = 1/2 for all the states relevant in this con¬ 
text. However, in atomic and molecular systems with 
more than one valence electron, S is at best an approx¬ 
imately good quantum number owing to inter-electronic 
interactions and spin-orbit coupling. Intercombination 
transitions in resonant energy transfer have been consid¬ 
ered in molecules [15] and quantum dots [16], and ex¬ 
perimental observations [13 HH] include important appli¬ 
cations in efficient organic light-emitting devices [S], In 
general, strong coupling to a bath of rotational, vibra¬ 
tional and motional states prevents the observation of 
coherent transport in these systems. 
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In this article, we examine the impact of singlet-triplet 
mixing on the long-range interaction between isolated 
ultracold atoms in the excitonic regime |19j . where de¬ 
phasing due to effects such as molecular motion is suffi¬ 
ciently small that the transport may be considered coher¬ 
ent. Our treatment may also find applications in molec¬ 
ular systems where coherent transfer plays a role. First 
we consider the effect of intercombination transitions on 
non-resonant van der Waals-type interactions, and sec¬ 
ond the effect of spin mixing on resonant hopping pro¬ 
cesses. We find that even in systems where spin mixing 
is relatively weak, near degeneracies between donor and 
acceptor states such as those that occur in Forster reso¬ 
nances can lead to novel transport effects which depend 
on the value of the spin quantum number S — here we 
refer to the actual electronic spin, rather than, e.g., to 
2-level excitation mapped to a pseudo-spin or to the spin 
magnetic quantum number Ms (be., the spin orienta¬ 
tion). Excitation transport in atomic lattices is an area 
of great current interest due to potential applications of 
cold-atom model systems to condensed matter research 
120. 21]. Here we consider a one-dimensional lattice of 
strontium atoms in 5 snd Rydberg states, thus adding 
strong inter-site interactions and extending the study of 
many-body systems beyond the two-level Hubbard model 
in cold atoms [22R24] . Ultra-cold Rydberg gases of diva¬ 
lent atoms are of growing interest in atomic physics [25r 
150] , and as well as systems where the precise details of the 
electronic wave function are known |3lti36] , they provide 
a route to precise control of the inter-particle spacing via 
optical lattices or tweezer arrays 

Since strontium has two valence electrons, two 5 snd 
Rydberg series with total angular momentum J = 2 
exist, one labelled as the singlet and one as a triplet. 
Early theoretical and experimental work showed that 
these energy eigenstates do not have a well defined spin 
due to their interaction with doubly excited “perturbers” 
of mixed singlet/triplet character 31, [32] 02] ■ These 
perturbers are coupled to the Rydberg states by inter- 
electronic interactions, resulting in a breakdown of LS 
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FIG. 1. (color online) (a) Two Rydberg atoms separated 
by a distance R are prepared in different spin states. Inter¬ 
actions between core and valence electrons (purple dashed 
arrow) leads to a breakdown of LS coupling, enabling an 
otherwise forbidden dipole-dipole interaction to take place 
(green dashed arrow), (b) The energy separation between 
the 130 3 D 2 , 30 3 D 2 ) pair states and the energetically closest 
spin allowed (red squares) and spin forbidden (blue circles) 
final pair states, showing a spin-forbidden resonance. 


coupling, which in turn affects the long-range inter¬ 
atomic interactions; the situation is depicted pictorially 
in Fig. [lja). We find that, for two interacting stron¬ 
tium atoms, the \n 3 D 2 ., n 3 D 2 ) states are close in en¬ 
ergy to the |(n — 2) ^ 3 , (n — 3) 3 F 3 ) states near n = 30 
(as shown in Figjljb)), which, combined with the spin¬ 
mixing in the Rydberg series, allows near-resonant trans¬ 
fer between two-atom states. The impact of these in¬ 
tercombination near-degeneracies on the van der Waals 
coefficients is considered in Section El For the case of 
four atoms, we find that a spin impurity (a 28 l F 3 atom) 
in a chain of 30 3 D 2 atoms can hop resonantly from site 
to site, showing a spin-forbidden propagation alo ng t he 
chain (as shown in Fig. [ 2 ] and discussed in Section III). 


II. MULTICHANNEL DIPOLE-DIPOLE 
INTERACTIONS 


We begin by considering two atoms in 5 snd 3 D 2 states. 
Theoretically, this situation is most easily treated using 
multi-channel quantum defect theory (MQDT), which 
provides a wave function for each single-atom energy 
eigenstate in terms of a superposition of LS'-coupled 
channels. Recently, we carried out an improved MQDT 
analysis of these states based on up-to-date experimental 
results, which gave the amount and nature of each elec¬ 
tronic state (singlet, triplet perturber) present in each of 
the J = 2 energy eigenstates [3T]. We use these wave 
functions to examine the long-range interaction between 
a pair of atoms prepared in the same 5 snd 3 D 2 energy 
eigenstate. Each of these pair states is coupled by electric 
dipole transitions to other final pair states (e.g. |P, P), 
|P, P), |P, F)). Because of spin mixing, the final states 
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FIG. 2. (color online) (a) Spin-forbidden propagation of 
a singlet spin impurity (28 Wa, initially located at Atom 1) 
along a short chain of triplet (30 3 D 2 ) Rydberg atoms, with a 
spacing of a = 2.0 pm between the atoms. Shading indicates 
the probability of each atom being in a S = 0 state as a 
function of time, (b) and (c) Time evolution of the probability 
of finding the atom in a state with spin quantum number 
S = 0 (b) and total angular momentum quantum number 
J = 3 (c) for different atoms along the four-atom chain. Line 
styles denote atom 1 (thick blue), atom 2 (medium green), 
atom 3 (thin black) and atom 4 (thin red). 
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may or may not differ in S from the initial state. A key 
parameter is the energy difference between the final state 
and the initial state — the so-called Forster defect - 
which must be compared to the strength of the coupling. 
Fig- 0b) shows an example where a near-degeneracy oc¬ 
curs in the spin-forbidden channel, i.e., where the spin 
labels of the initial and final states are different. Thus, al¬ 
though the spin mixing, and hence the coupling, is weak, 
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this spin-forbidden process can become important. 

More concretely, to describe the long-range interac¬ 
tions, we consider each atomic energy eigenstate state, 
T, to be a sum over the MQDT channel states, (j>k, such 
that A-k4>kXk (where Xk is a function describing 

the angular, spin and remnant core state wave functions 
eh). Using these state vectors, the long-range inter¬ 
actions can be calculated either perturbatively or non- 
perturbatively (by diagonalizing an effective Hamiltonian 
matrix in a basis of pair states [IS]). The values of the 
coefficients A/-, as well as numerical dipole matrix ele¬ 
ments, are provided in }31j . Throughout this paper, we 
only consider atoms that are initially in a stretched state 
(J = |Mj|), with the internuclear axes of the interact¬ 
ing atoms being aligned with the z-axis. Stretched states 
do not have any degeneracies in Mj 1 + Mj 2 , thereby re¬ 
ducing the number of states that need to be considered 
(even allowing for the fact that the dipole-dipole interac¬ 
tion couples stretched states to non-stretched states). 

The dipole-dipole interaction Hamiltonian for two 
atoms with the internuclear axis aligned with the z-axis 
is given by 03] 04] 

A 'jr 

H dd = ~YR? Tir2 ( y bi(^i) F i,-i( f 2) 

+ Yi_i(fi)yi,i(f 2 ) + 2Fi ! o(ri)yi,o(f ; 2)), 


where R is the interatomic distance, r\ and 7' 2 are the 
radial electronic coordinates for atoms 1 and 2, rq and 
f 2 the angular electronic coordinates and Y[ rn denotes a 
spherical harmonic. The matrix elements of 0 are thus 
products of angular factors and dipole matrix elements 
for each atom, with a R~ 3 dependence. 

Treating the dipole-dipole interaction using perturba¬ 
tion theory the first-order expression vanishes for two 
atoms in the same state. Second-order perturbation the¬ 
ory results in a C e R~ 6 interaction, with [13] 
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where atom 1 is in state and atom 2 is in state 

l’^ 2 °' ) ) ( n °t e that these are not necessarily the same 
states), and the sum over i runs over all the pair states 
dipole-coupled to By using the MQDT expan¬ 

sions of the wave functions in terms of channels and us¬ 
ing the channel fractions and doubly excited state dipole 
matrix elements from EH, Eq. ([2]) can be evaluated nu¬ 
merically. 

The resulting values of the Cq coefficients for a pair 
of Sr atoms both in the same 5 snd 1 _D 2 ( \Mj\ = 2) 
state or the same 5 snd 3 D 2 ( \Mj\ = 2) state are 
shown in Table [I] The contributions from “spin-allowed” 
(i.e. singlet-singlet, triplet-triplet) and “spin-forbidden” 
(singlet-triplet) intermediate pair states are presented 
in Fig. [3] Large contributions from singlet-triplet pair 
states are found in both series around n = 16 where the 


TABLE I. The Cq coefficients for the singlet and triplet 
J = 2 5 snd (|A/j| = 2) configurations of Sr (where n denotes 
the principal quantum number). The coefficients are given 
in atomic units, where C6(GHz / 1 m 6 ) = 1.4448 x 1CF 19 C6 
(atomic units). 
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effect of the 4d6s 1 D 2 and 3 D 2 perturbers is at its max¬ 
imum [321 03 ] ■ The overall Cq coefficients for states in 
this region differ significantly from predictions based on 
single-channel quantum defect calculations for Rydberg 
states below n = 30 05]. Above n = 30, however, the 
calculated single-channel and multichannel values differ 
by less than 2% of the overall Cq, thus validating the use 
of a one-electron treatment for high-lying Rydberg states 
of strontium 3Ti 03 ] 06 ]. 

Also visible in Fig. [3] is a large singlet-triplet con¬ 
tribution for 3 D 2 states close to n = 30. This arises 
due to the Forster resonance in the |n 3 D 2 ,n 3 D 2 ) —» 
|(n — 2) 1 J 7 3, (n — 3) 3 F 3 ) channel shown in Fig. [ljb). The 
uncertainties in the energy levels used to calculate the 
Cq coefficients 03] are large enough that the location 
of the Forster resonance can change by one value of n; 
however the resonance is always present to within the 
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FIG. 3. (color online) Absolute value of the contributions 
to the C 6 coefficients in atomic units (with the dominant n 11 
scaling factored out for clarity) from the singlet-singlet (light 
blue diagonally-lined bars), singlet-triplet (dark green) and 
triplet-triplet (light green horizontally-lined bars) pair states 
acting as intermediate states. The initial states are taken to 
be in their stretched state, with J = Mj. 


error of these energy level measurements [32] 03, I TT] . 
The small Forster defect in this channel means that 
second-order perturbation theory breaks down at rela¬ 
tively large interatomic distances. We therefore turn to 
a non-perturbative calculation. Fig. [4] shows the non- 
perturbative Born-Oppenheimer potential curves in the 
vicinity of the |30 3 Z?2, 30 3 D 2 ) asymptote, which has a 
spin-forbidden avoided crossing at relatively large dis¬ 
tances (R ~ 0.5 /jin) with the 128 X F 3 , 27 3 F 3 ) asymptotic 
pair state. Without the mixing between the triplet and 
singlet series, the avoided crossing would not exist. While 
this Forster defect (522 MHz) is not small compared to 
that found in alkali atoms [35H5U] . the key point here is 
that it is much smaller than the defect for the dipole- 
allowed pair states (8.84 GHz). As a result, the inter¬ 
action between the spin-forbidden pair states is stronger 
than could be expected in view of the smallness of the 
singlet-triplet mixing in these Rydberg states. 


III. SPIN CHAIN OF STRONTIUM RYDBERG 
ATOMS 

As another illustration of the impact of this intercom¬ 
bination Forster resonance on resonant energy transfer, 


we now examine the propagation of a singlet “impurity” 
in a short chain of four equally spaced atoms. While 
studies of state transport in lattices of Rydberg atoms 
have already been carried out ED Eg, as far as we know 
the impact of singlet-triplet mixing in this context has 
not been previously considered. Denoting the 30 3 D 2 , 
28 1 F 3 , 27 3 F 3 , and 28 3 F 3 states by |0), |1), |2), and |3), 
respectively, we numerically calculate the time evolution 
of the system at time t after the 11000) state is pre¬ 
pared. Restricting the dynamics of each atom to these 
four states is justified by the fact that the C 6 coefficient of 
the |30 3 D 2 , 30 3 D 2 ) is dominated by the Forster-resonant 
|28 3 F 3 , 27 3 F 3 ) and the non-resonant |28 3 F 3 , 27 3 F 3 ) pair 
states. All other pair states contribute less than 15% to 
the C 6 coefficient of the |30 3 D 2 , 30 3 D 2 ) state and are far 
enough away in energy to be neglected. Fig. [4] shows the 
consequence of only choosing the four single atom states, 
30 3 Z? 2 , 28 3 F 3 . 27 3 F 3 . and 28 3 F 3 . The potential curve for 
the 1 30 3 D 2 , 30 3 D 2 ) asymptote is well reproduced. We in¬ 
clude all values of Mj that contribute. 

In order to perform the time-dependent calculation of 
the spin chain, we write the total Hamiltonian as H = 
iJo + ifint, where Ho is the Hamiltonian of the four atoms 
with infinite lattice spacing. The interaction Hamiltonian 
Hint can be written as 


Hint = 




where i and j label the many-body states and the matrix 
element V is given by 

v ji =X]^l i) ^2 ) ^3 ) ^4 ) |flpql^l j) ^2 j) ^3 i) ^4 ) )- ( 4 ) 

p<q 


Here .ff^^is the dipole-dipole interaction Hamiltonian 
given in pi) between lattice sites p and q. 

As the dipole-dipole interaction is a time-independent 
perturbation, we expand the eigenstates of H, 
|$ a (o,t)), in terms of the four-atom states of the 
lattice with infinite spacing, such that |4> Q (a, t)) = 
exp(ie a (a)t)^-L^ ) (a)|^^ ) ^ ) ^ ) ^ ) ), where e a (a) 
are the eigenenergies of H. The initial state vector 
I’F^'F^ 0 ' 'k 3 0 ' 4, 4 °^) is projected from the basis of the bare 
pair states into the eigenbasis. The exp(ze a f) factors are 
then easily determined, and the final state vectors are 
projected back into the original basis. To calculate the 
probabilities of the spin chain to be in a state i, we use 
the square magnitude of the coefficients 


c W (a, t) = Y ( U a ] ) exp (ie a t) . (5) 

a. 

Fig- m shows the evolution of the probability of each 
atom being in a singlet state for a lattice spacing of 2 /jrn 
(a spacing that can be engineered using two crossed 1550 
nm laser beams [231153]). The spin can be seen to propa¬ 
gate along the chain of atoms and back, although there is 
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FIG. 4. (color online) Dipole-dipole potential curves for two 
atoms around the state labelled |30 3 Z)2, 30 3 I?2) at infinity. 
The blue dashed curves represent the full non-perturbative 
calculation including all the relevant pair states, and the solid 
curves show the results of the four-level approximation (the 
red curve highlights the initial state, 130 3 Zi> 2 , 30 3 I?2)). Only 
the Mj x = Mj 2 = 2 states are shown for the |30 3 T>2, 30 3 D2) 
state. 


additional state transfer due to competing second-order 
interactions. The calculation shown in Fig. [2] includes 
the interactions between all the atoms, not just nearest- 
neighbour interactions. Nevertheless, a clear propagation 
of a spin singlet state through the chain can be seen, a 
phenomenon that can only occur due to spin-mixing. 

For the parameters of Fig. [2j the dynamics arise pri¬ 
marily from the spin-forbidden dipole-dipole coupling be¬ 
tween the four linearly independent states of the 0, 0, 0,1 
family, namely the states which reduce to linear combi¬ 
nations of the 11000), |0100), |0010) and |0001) states in 
the limit of infinite lattice spacing. As states 2 and 3 are 
unimportant here, there is essentially no difference be¬ 
tween the probability for a particular atom of the chain 
being in a S = 0 state, P(S = 0), and the probability 
of it being in a J = 3 state, P{J = 3) (compare parts 
(b) and (c) of Fig. [2]). However, this is not the case 
for smaller lattice spacings. As shown in Fig. [5ja) , the 
four 0,0, 0,1 states exhibit avoided crossings between 1.2 
and 1.6 /mi with the 0,1,1,2 family of states, namely 
states in which, for a —> oo, one of the four atoms is in 
state 0, two are in state 1 and one is in state 2. Due to 
these spin-forbidden Forster resonances and to the larger 
strength of the dipole-dipole interaction, the dynamics 
of the chain at a = 1.35 fim is more complex than at 
a = 2 /mi [Figs. [5](b) and (c)]. In particular, for some 
of the atoms the stronger coupling with state 2 at this 
smaller lattice spacing results in striking differences be¬ 
tween P(S = 0) and P(J = 3). The key feature remains 
that the spin-orbit and inter-electronic effects responsi¬ 
ble for the breakdown of LS coupling within each atom 
manifest spatially in the collective state of the spin chain. 

In regards to an experimental study of this intercom¬ 
bination dynamics, we note that the natural lifetimes 



a (//m) 



0 20 40 60 80 100 


Time (ns) 



0 20 40 60 80 100 


Time (ns) 

FIG. 5. (color online) (a) Born-Oppenheimer energy curves 
for a spin chain of four atoms equally separated by a distance 
a from their nearest neighbours, with the energy of one of the 
states of the 0, 0,0,1 family (solid black curve) and one of the 
states of the 0,1,1, 2 family (dashed black curve) highlighted 
for clarity. Subfigures (b) and (c) show the probability of 
atoms in the spin chain to be in a S = 0 (b) and a J = 3 
state (c), for a lattice spacing of a = 1.35 //m. Line styles 
denote atom 1 (thick blue), atom 2 (medium green), atom 3 
(thin black) and atom 4 (thin red). 


of the 30 3 Z?2 and 28 ^ states (2.3 /is and 6.9 /is, re¬ 
spectively eh m) are much longer than its time scale. 
The triplet F states lifetimes are unknown, but can 
be expected to be similar in magnitude to that of the 
30 3 D -2 and 28 ^F :i states. A possible experimental reali¬ 
sation could employ a strontium quantum gas microscope 
(QGM) [55] [56], where atoms are loaded into the lowest 
vibrational band of a 3-D lattice using a Mott-Insulator 
transition EH, and a high numerical aperture objective 
provides the required single-site readout. Single-site de¬ 
tection of Rydberg atoms in a QGM has already been 
demonstrated [1581 . Since the atoms are in the lowest 
band, dephasing effects due to the uncertainty in the ini¬ 
tial positions @001531 are minimised. Localisation to 
< 50 nrn is possible for reasonable lattice depths ,55] . 
which is less than the width of the avoided crossing in 
Fig-0 and which is sufficient to observe coherent trans- 
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port under the conditions of Fig.[2j We note that due to 
the fast timescales, the lattice could be switched off dur¬ 
ing the transport process. Finally, in order to image the 
state of the chain, short microwave pulses could be used 
to state-selectively transfer the population to other Ry¬ 
dberg states that do not interact resonantly, thus “freez¬ 
ing” the dynamics. 

IV. CONCLUSIONS 

In conclusion, we have shown that intercombination 
transitions in Sr Rydberg atoms not only lead to a break¬ 
down of LS coupling but also allow dipole-forbidden exci¬ 
tation hopping along a chain of atoms via resonant long- 
range dipole-dipole interactions. We find that intercom¬ 
bination Forster resonances can have a substantial im¬ 


pact on long-range interaction. They can also lead to 
spatially separated dynamics between spin angular mo¬ 
mentum and total angular momentum. Although we use 
Sr Rydberg states as an example, the ubiquity of spin 
mixing makes it likely that other systems may also show 
similar effects. 
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